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We determine the scaling functions describing the crossover from Ising-like critical behavior to 
classical critical behavior in two-dimensional systems with a variable interaction range. Since this 
crossover spans several decades in the reduced temperature as well as in the finite-size crossover 
variable, it has up to now largely evaded a satisfactory numerical determination. Using a new 
Monte Carlo method, we could obtain accurate results for sufficiently large interactions ranges. Our 
data cover the full crossover region both above and below the critical temperature and support 
the hypothesis that the crossover functions are universal. Also the so-called effective exponents are 
discussed and we show that these can vary nonmonotonically in the crossover region. 
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I. INTRODUCTION 



The crossover from Ising-like to classical critical behavior has attracted renewed attention in recent years. This 
crossover behavior occurs in many thermodynamic systems, such as ionic solutions, simple fluids, fluid mixtures, and 
polymer mixtures. The Ginzburg criterion fjl] states that sufficiently close to the critical point these systems exhibit 
critical exponents belonging to the 3D Ising universality class. At larger distances from the critical point, but still 
within the critical region, classical (mean-field-like) critical exponents are observed. Although this appears to be 
a well-established picture, the precise nature of the crossover between these two universality classes is still subject 
to investigation. For example, Anisimov et al. recently claimed || to have observed an "effective" susceptibility 
exponent which varied nonmonotonically from its classical value 7mf = 1 to its Ising value 71 ~ 1.24 when the critical 
point was approached. Later, the possibility of such behavior within the critical domain was questioned by Bagnuls 
and Bervillier, see Refs. HQ. On the other hand, Fisher has argued || that nonmonotonical variation of effective 
critical exponents is not necessarily an indication of nonuniversal behavior. Other questions concern the size of the 
crossover region, which is expected to span several decades in the crossover variable j(| , and the size of the temperature 
region around T c within which Ising-like behavior is observed [Q. Until now it has turned out to be very difficult to 
accurately observe the full crossover region in numerical simulations. A major effort has been undertaken in Ref. j^] 
for three-dimensional polymer mixtures, where crossover occurs as a function of the polymer chain length. However, 
despite chain lengths of up to 512 monomers, the results did not span the full crossover region. For this reason, 
Mon and Binder J^] turned their attention to the two-dimensional Ising model with an extended range of interaction, 
where a crossover from Ising-like to classical critical behavior occurs when the range R of the spin-spin interactions is 
increased (suppressing the critical fluctuations). In two dimensions one can not only access larger interaction ranges, 
but also both asymptotic regimes are known exactly and the variation of the critical exponents is considerably larger 
than in the crossover from 3D Ising-like critical behavior to classical critical behavior. Mon and Binder derived the 
(singular) R dependence of the critical amplitudes of scaling functions and carried out Monte Carlo simulations to 
verify these predictions numerically. Even in these two-dimensional systems, the mean-field regime turned out to be 
only barely reachable. 

In a recent paper p0[ , we rederived the predictions of Mon and Binder from renormalization theory and also 
obtained the R dependence of various corrections to scaling, such as the shift of the critical temperature with respect 
to the mean-field critical temperature. Furthermore, larger interaction ranges and system sizes were accessible to our 
numerical simulations thanks to a dedicated Monte Carlo algorithm. This enabled us to actually verify the theoretical 
predictions in two-dimensional systems. In this paper, we show that the simulations presented in Ref. |l(]| ] allow a 
full mapping of the finite-size crossover curves for various quantities. However, these curves describe the finite-size 
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dependences of critical amplitudes, which have (to our knowledge) not been observed experimentally. Therefore 
we have also carried out simulations at temperatures further from the critical temperature in order to observe the 
thermal crossover of these quantities. The results of these simulations, which partially have been reported in an earlier 
note are presented as well. The fact that in our model both the temperature distance from the critical point and 
the interaction range can be varied turns out to be essential to observe the full crossover region. 

The outline of the remainder this paper is as follows. After a short recapitulation of the model under investigation 
(Sec. |J) we start in Sec. Ill with finite-size crossover scaling. We discuss the required system sizes and interaction 
ranges and obtain the crossover curves for the absolute magnetization density, magnetic susceptibility, and the spin- 
spin correlation function over half the system size. Thermal crossover scaling is treated in Sec. IV, where we consider 
the approach of T c both in the symmetric phase (T > T c ) and in the state of broken symmetry (T < T c ). Again, 
crossover curves are obtained for the order parameter and the susceptibility. The various aspects of these curves are 
discussed in some detail. Graphs of the logarithmic derivatives of the crossover curves, which can be associated with 
so-called effective critical exponents as measured in experiments, are presented in Sec. [y|. In Sec. VI we end with a 
summary of our conclusions. 



II. SHORT DESCRIPTION OF THE MODEL 

Let us first briefly recall the model as it was introduced in Ref. |^| . This is a two-dimensional Ising system consisting 
of L x L lattice sites with periodic boundary conditions. Each spin in the system interacts equally with its z neighbors 
lying within a distance R m . This defines the coupling between two spins Sj and Sj at a distance r as 

K ij =K(r) = (f™ ?/zv n (!) 



l] K '-\0 if r>R 
In the absence of an external magnetic field the Hamiltonian is 

H/fc B T=-^^X(|r l -r J |) SlSj , (2) 



i j>i 



where the sums run over all spins in the system and denotes the position of spin Si . To suppress lattice effects we 
use an effective interaction range R, defined as 

2 _ r,fl<jj 



Rf = 



- ^ \T t - Tj | 2 with \n -rj| < Rm . (3) 



z . 

3^ 



For large ranges, R approaches the limiting value R m /y/2. 



III. FINITE-SIZE CROSSOVER SCALING 



A. General considerations 

It has been shown by Binder and Deutsch C3] that crossover scaling can be combined with finite-size scaling by 
including the dependence on the crossover variable in the probability distribution function of the order parameter. 
Indeed, just as crossover in the thermodynamic limit is described as a function of the reduced temperature divided by 
the Ginzburg number, it can be described as the function of a size-dependent crossover variable G in finite systems. 
In Ref. P], this crossover variable was derived as G = LR~ A t^~ d \ where L is the linear system size and d denotes 
the dimensionality. This also follows from the renormalization treatment in Ref. [ jTo| . In short, sufficiently close to 
the Gaussian fixed point (i.e., for a sufficiently large interaction range R) the critical behavior will be classical. In 
terms of a renormalized Landau-Ginzburg- Wilson (LGW) Hamiltonian in momentum space, this implies that the 
coefficient of the </> 4 term must be much smaller than that of the 2 term, uL 4 ~ d /R 4 <C 1 (cf. in particular Eq. (6) 
of Ref. @), which again leads to the crossover parameter LR 4 / ( 4 ^ , where for the moment we assume that u is of 
order unity. 
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In Ref. (h}, we focused our attention on the critical finite-size amplitudes in the limit of L — > oo. Here we will 
examine the crossover in the corresponding data for finite system sizes. Since the crossover regime is expected [^[(^] 
to span several decades in the crossover variable G = L/R 2 , it is numerically not feasible to observe both asymptotic 
regimes by merely varying the system size L while keeping the range R fixed. Therefore we construct the curve by 
combining the results for various values of R, cf. Ref. a]. Indeed, the Ising regime (L/R 2 ^> 1) is easily reachable, 
although the results for very small ranges do not conform well to the leading R dependence of the critical scaling 
functions Q and are thus, at first sight, not well suited for constructing the crossover curve. The mean-field regime 
(L/R 2 <C 1), however, poses more substantial problems. If the linear system size L is made too small, the numerical 
results exhibit strong finite-size effects. Therefore L must be at least of the order of the interaction range. More 
precisely, boundary effects will occur for systems for which L w R m and the smallest possible value of the crossover 
variable G is roughly equal to R m /R 2 y/2/R. Thus, large ranges are required to reach the regime where G< 1. 
In a conventional Monte Carlo algorithm, the efficiency of simulations rapidly decreases with increasing interaction 
range. This limitation has been circumvented by applying a dedicated cluster algorithm, as explained in Ref. p0[ . 
Still, a problem remains. Namely, the finite-size crossover scaling is valid at the critical temperature. Any deviation 
from this temperature will lead to systematic errors in the analysis. Since the (range-dependent) critical temperatures 
are determined in the Ising limit, i.e. from system sizes L > R 2 , large interaction ranges require very large system 
sizes for an accurate determination of T c . For example, the most efficient way to obtain data for G ~ 0.02 is to 
simulate a system with L = 100 and R m = 100 (R w 70). However, an accurate determination of T C (R = 70) requires 
system sizes of at least L — 5000, whereas we have carried out simulations for system sizes up to 1000 x 1000 lattice 
sites. This has been solved as follows. The renormalization treatment in Ref. [|l0| predicts the form of the function 
describing how T C (R) deviates from the mean- field critical temperature when R varies. By fitting this function to the 
accurately determined critical temperatures in our previous study an expression is obtained for T C (R) from which the 
critical temperatures for very large ranges can be calculated to a relatively high accuracy. The shift of T c is expressed 
by 

T c =T c MF + ^[l + a 2 lni? 2 ]+||, (4) 

where T C MF = 1 and the last term is a higher-order correction omitted in Ref. p0[ . A least-squares fit for 16 < R 2 < 70 
(32 < R 2 m < 140 in Ref. @) yielded a x = -0.267(6), a 2 = 1.14(3), and a 3 = -0.27(3). Figure | shows the critical 
temperatures and expression (0) with the appropriate coefficients. 



B. Absolute magnetization density 

In the Ising regime, the absolute magnetization density scales (at criticality) asymptotically as (\m\) = L~ 1 / 8 do(R), 

where the critical amplitude do is a function of R, do oc i?~ 3 / 4 . In the mean-field regime (\m\) does not depend on R, 

but is simply proportional to L~ x l 2 . When plotting (\m\) as a function of G = L/R 2 a data collapse is obtained if 

it is multiplied by a factor L x i? - ^ -1 ). This resulting quantity is proportional to G x ~ 1 ^ s in the Ising regime and to 
qx-i/2 in the 

mean-field regime. A suitable choice is x = 1/2, because this yields a quantity which is still independent 
of R in the mean-field regime. Indeed, it is shown in the Appendix that in a two-dimensional system in which all 
spin-spin interactions are equally strong, 

and (|m|)yT will thus approach 12 1 / 4 r(l/2)/r(l/4) = 0.909890588... in the limit of G -> 0. Remark that our 
requirement L > y/2R unambiguously relates the limit G — > to the mean- field (R — > oo) limit. In Fig. |2|(a) we 
have plotted the absolute magnetization density multiplied by the square root of the system size versus the crossover 
variable. Interaction ranges from R 2 n = 2 to R 2 n = 10000 were included, where the data for R 2 ^ = 5000 and 
R^ = 10000 (spanning the range 0.02 < G < 0.2) have been obtained at temperatures calculated from Eq. (Q): 
K c (R m = v^OOO) = 6.3746(3) x 10" 5 and K c (R m = V10000) = 3.18491(9) x 10~ 5 . The crossover curve evidently 
spans approximately three decades in G. In the limit of G — ► it gradually approaches a horizontal line. For 
G> 1 the picture is not very clear. The data points for each single value of R lie on a straight line with slope 3/8, 
corresponding to the Ising asymptote, but the asymptotes only coincide for large ranges (cf. Fig. 4 in Ref. ||). The 
reason for this is that, as mentioned above, for small ranges the critical amplitudes do not conform to the leading 
i?~ 3 / 4 dependence. This can be cured by invoking the renormalization treatment of Ref. |jTc| . Indeed, the theory 
predicts the structure of the corrections to the leading R dependence of the critical amplitude, 
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d = b R~ 3/4 



R 2 



(b 1 +b 2 lni? 2 



(6) 



This "finite-range correction" is very similar to the shift of the critical temperature in Eq. (|4|) but originates from a 
different term in the renormalized LGW Hamiltonian. To illustrate this correction graphically, we have reproduced 
Fig. 6 from Ref. |hJ] and included the result of a least-squares fit of Eq. @ to the data, see Fig. |. The curve clearly 
yields an excellent description of the critical amplitudes, even for small ranges. We have used this fit to construct 
a clear crossover curve for the magnetization density on which the data for all values of R collapse. To this end, 
all data are divided by the correction factor between square brackets in Eq. (^). The result is shown in Fig. |^(b). 
One observes that in the Ising regime all data perfectly collapse on a common asymptote with slope 3/8. For G 
small the data indeed approach the mean- field prediction (|). The fact that at G « 0.2 the data for R 2 m = 5000 and 
R^ = 10000 coincide with those for i?^ = 72, 100, 140 confirms that the critical temperatures for the large ranges 
have been estimated accurately. The center of the crossover region lies between G = 0.1 and G = 1.0 which shows 
that the parameter u is indeed of order unity. Finally, it is particularly encouraging that no remaining finite-size 
effects, causing deviations from the curve, are visible in Fig. |2|(b), despite the factor that the correction factor was 
calculated in the L — > oo limit and hence does not compensate for such higher-order finite-size effects. 



C. Magnetic susceptibility 



The procedure described above for the absolute magnetization density can be applied to the magnetic susceptibil- 
ity x, which we have calculated from the average square magnetization, \ — L d (m 2 ). At T — T c , the susceptibility is 
in the Ising regime proportional to L 7 / 4 i?~ 3 / 2 , and in the mean- field regime it scales proportional to L. To obtain a 
data collapse for x as a function of G, one has to multiply the finite-size data by L x R~( 2x+2 \ where a suitable choice 
is given by x = — 1. In the mean- field limit, \ jL approaches VT2r(3/4)/r(l/4) = 1.17082866 . . . (see Appendix). As 
shown in Ref. fl0|| , the deviation from the leading range dependence of the critical amplitude is very similar to that 
of the absolute magnetization density, 



Pa = <7o 



R -3/2 



(7) 



where now one additional higher-order correction is required. Therefore we only show the resulting crossover curve 
for the susceptibility after the data have been divided by the correction factor between square brackets, see Fig. ||. 
Again, both the mean- field asymptotic result and the Ising asymptote (slope 3/4) are clearly reproduced, with a 
perfect collapse for all ranges. 



D. Spin— spin correlation function 



Closely related to the magnetic susceptibility is the spin-spin correlation correlation function g(r). In our simula- 
tions we have sampled g(L/2), which scales both in the Ising regime and in the mean-field regime as x/L 2 ■ Thus, 
we obtain a data collapse by multiplying the finite-size data by L x R~( 2x ~ 2 \ in which we have set x = 1. After 
correcting for the higher-order range-dependent corrections in the critical amplitude [which have the same structure 
as those in Eq. (0)] we obtain the graph shown in Fig. The full crossover curve can be mapped and shows a close 
resemblance to that for the susceptibility, including the approach of the asymptotic mean-field value. In the range 
0.2 < L/R 2 < 1.0, the data do not precisely coincide on a smooth curve. This is due to nonlinear finite-size effects, 
which are for the spin-spin correlation function apparently larger than for the absolute magnetization density or the 
magnetic susceptibility. We will pay more attention to these deviations when discussing the universal amplitude ratio 
(see below). It should be noted that the critical amplitudes listed in Table V of Ref. (lO) have to be multiplied by a 
factor 2- 1 / 4 = 0.84089642 . . . in order to obtain the correct values. 



E. Universal amplitude ratio 

The amplitude ratio Ql = (wL) 2 /( m i) * s a size-dependent quantity, which takes a universal value Q in the L — > oo 
limit. That is, it is calculated by taking the ratio of the square of the magnetization density and the fourth power of it 
in a finite geometry and subsequently taking the limit L — > oo. For T > T c , Q approaches the Gaussian value Q = 1/3 
and for T < T c it approaches the maximum value Q = 1. At criticality, the amplitude ratio is known exactly 
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in the mean-field case, Qmf = 0.45694658 . . . HQ and to a high accuracy in the two-dimensional Ising model, 
Qi « 0.856216(1) In Ref. Jioj, Ql(K c ) was plotted for a large interaction range (R 2 n — 140) as a function of the 
system size. The approach of the Ising value was clearly visible for L large, but for small system sizes Q first decreased 
towards Qmf and then started to show strong finite-size effects. Evidently, it is a better approach to construct the 
true crossover curve for Q(K C ) by plotting finite-size data for Q for various ranges versus the crossover variable. This 
is shown in Fig. ^(a). Several remarks apply to this graph. Firstly one notes that L/R 2 is indeed the appropriate 
crossover variable: a reasonable collapse is obtained for all values of L and R. However, some remarkable deviations 
from this scaling behavior are present, which are most clearly visible in the range 0.2 < L/R 2 < 0.6, but also present 
around L/R 2 = 10. Similar effects were already observed in the spin-spin correlation function, but now the effects 
stand out much more pronounced because we have employed for the amplitude ratio a linear instead of a logarithmic 
vertical scale. These deviations are due to nonlinear finite-size corrections, as can be seen clearly by zooming in 
into the deviations, see Fig. 0. The data points for R 2 m = 5000 and R 2 m = 10000 may serve as a reference for the 
location of the "true" crossover curve. One observes that for each of the ranges i? 2 ^ = 72, 100, and 140 the deviations 
from this curve increase with decreasing system size, which indeed shows that the effects are caused by finite-size 
corrections. If the deviations had been caused by e.g., an inaccurate determination of the critical temperature, the 
effects would have increased with increasing system size. Unfortunately, it is not easy to separate these corrections 
from the leading crossover behavior (except graphically), unless the full crossover function is known (which in turn 
would limit the use of a numerical determination). Of course this problem can be circumvented by determining the 
crossover at these values for G from systems with a larger system size and a larger interaction range. The deviations 
around L/R 2 = 10 are caused by the same effect, but now for systems with small R. Although the amplitude ratio 
is more sensitive — even if one takes into account the difference in scale — to these finite-size effects than (m 2 ) = x/L 2 
and (to 4 ) individually (the curve for the latter is not shown here, but its smoothness is comparable to that of the 
susceptibility), Q is less sensitive to corrections to the leading range dependence. Indeed, for (to 4 ) these corrections 
are again of the form [1 + R~ 2 (si + s 2 lni? 2 ) + R^ 4 s 3 ] and Q must thus be divided by 

[1 + R- 2 ( qi + g 2 \nR 2 ) + R-%} 2 

l + Rr 2 ( Sl + s 2 h-iR 2 ) + R- 4 s 3 • 1 ' 

The coefficients si, s 2 , and S3 have been determined from a least-squares fit to the critical amplitudes of (m 4 ) and 
qi, q 2 , and q 3 come from Eq. (0). Figure || shows the correction factors for (m 2 ), (to 4 ), and Q. Evidently the latter 
factor (H) is much closer to unity than the former two. Figure ||(b) shows Ql(K c ) divided by the correction factor (||), 
which indeed shows only slightly less scatter than the graph without this correction factor. In particular the deviations 
for the larger ranges do not disappear. 

IV. THERMAL CROSSOVER SCALING 
A. General considerations 

The finite-size crossover scaling studied in the previous section is an intrinsic finite-size effect which is not observable 
in thermodynamic systems. For this reason it is important to study its temperature-dependent counterpart as well. 
This so-called thermal crossover, which was from a phcnomcnological scaling point of view already considered in 
Ref. (l(|, is of course closely related to finite-size crossover: in finite systems crossover to mean-field-like behavior 
occurs when the system size has been decreased to the appropriate power of the interaction range (i.e. L ~ 
or L ~ R 2 for d = 2), whereas in the temperature-dependent case this crossover occurs when the temperature distance 
to the critical point is such the correlation length has become of the order of an appropriate power of the interaction 
range. In the latter case, the precise crossover location is determined by the Ginzburg criterion, t^^^R^vT 1 sa 1, 
where u is the coefficient of the 4> A term in the LGW Hamiltonian. It should be kept in mind that these considerations 
are valid only within the critical region, i.e. care must be exercised to keep the reduced temperature sufficiently small. 
When studying thermal crossover in practical simulations one has the additional complication that sufficiently close 
to T c the correlation length will always be bounded by the finite system size, which is precisely the situation one 
wants to avoid. So relatively large system sizes are required. 

As follows from the Ginzburg criterion, the appropriate scaling variable in two dimensions is tR 2 and one can 
therefore study thermal crossover effects by varying the interaction range as well. This is essential because of the 
following. For small values of R, t has to be made rather large to cross over to classical critical behavior and it is 
possible that one leaves the critical region before reaching the classical regime. On the other hand, if one only studies 
systems with large interaction ranges, t has to be made very small to observe Ising-like critical behavior. However, 
for such small values of t extremely large system sizes are required to avoid finite-size effects. Therefore we have 
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constructed, just as in the previous section, crossover curves from results for various ranges. We have carried out 
simulations for the interaction ranges studied in Ref. ||ic|| at temperatures further below T c and also generated data 
for the interaction ranges R 2 ^ = 500, 1000, 4000, and 10000. Table | summarizes some properties of these systems. 
Simulations have been carried out down to temperatures as low as T 0.5T C . For the order parameter crossover can 
only be studied in the phase of broken symmetry, but for the susceptibility we have also considered the symmetric 
(T > T c ) phase. Since in this phase no saturation effects occur, much smaller interaction ranges suffice to span the 
full crossover region, as we will show below. 



B. Absolute magnetization density 

As derived in Refs. P,M the absolute magnetization density scales, sufficiently close to the critical point, as 
(\m\) oc (— t}P R( 2d P- d )/<J-^f (t < 0), which for the two-dimensional case yields (\m\) oc (— t) 1 / 8 ^ -3 / 4 . In the mean- 
field regime, on the other hand, the magnetization density is simply proportional to (— t) 1 / 2 . When plotted as a 
function of tR 2 , a data collapse for all ranges is now obtained if the magnetization density is multiplied by R. 
Figure ||(a) shows the corresponding plot. We will discuss the various aspects of this graph in some more detail. The 
overall picture suggests that the data roughly follow the Ising asymptote (slope 1/8) for small values oitR 2 and then 
gradually approach the mean-field asymptote (slope 1/2) for large values of tR 2 . Here "small" and "large" refer to 
the absolute value of tR 2 and "slope" is generally used for the logarithmic derivative, dln(|m|)/dln \t\. For very small 
values of tR 2 the data start to deviate from the Ising asymptote at an L-dependent location and approximately follow 
(for temperatures closer to T c ) a horizontal line. Here one has entered the finite-size regime, where the correlation 
length is limited by the system size. This is the case which was studied in the previous section. The width of 
this regime depends (for general d) both on the system size and the interaction range, as can be read off from 
the universal scaling functions derived in Ref. ficfl . Indeed, the temperature-dependent argument of these functions 
is tL Vt R~ 2 ( 2vt ~ d ^( 4L ~ d } (y t = 1 in the 2D Ising universality class) and the width of the finite-size regime is thus 
proportional to L~ Vt R 2 ^ Vt ~ d >' = L _1 . Note that the absence of any range dependence is not a general feature 
and even for the two-dimensional Ising model only true to leading order (cf . Fig. 5 of Ref. [jl0| ) . Higher order terms 
will entail range-dependent factors that involve (for d = 2) logarithms of R. Outside of the finite-size regime, the data 
for each individual range first lie approximately on the Ising asymptote, which has been drawn with an amplitude 
such that it coincides with the data for R 2 m = 2. For the smaller ranges the amplitudes of the asymptotes show 
a considerable range dependence, whereas for larger ranges the amplitudes converge. Upon further decrease of the 
temperature (increase of the absolute value of t) several types of behavior occur: for the smallest range (R 2 n — 2) the 
data points still lie on the Ising asymptote. For R 2 n = 4 and R^ = 10 the data leave the Ising asymptote at sufficiently 
low temperatures and then follow a nearly straight line with a slope that lies between the Ising and the mean-field 
asymptote. In these cases one has left the critical region without ever reaching the asymptotic mean-field regime. 
For each range the data for all system sizes coincide, as they should outside of the finite-size regime. For = 72 
and R 2 m = 140 the mean-field asymptote is approached much closer. However, if the temperature is decreased further 
below the critical temperature the data points start to deviate from the asymptote again. This effect is caused by 
saturation of the magnetization and can be quantitatively described with mean-field theory, as we will show below. 
Turning to even larger ranges, we see that the data now really reach the asymptote with slope 1/2 and follow it for 
up to one decade in the crossover variable (for the largest range we have studied) before saturation sets in. Also 
the exact amplitude -\/3 (see below) of the asymptote is precisely reproduced, which shows again that the critical 
temperatures of the systems with large interaction ranges have been accurately determined: a deviation would have 
shifted the graph along the horizontal axis. 

We will now first consider the offset of the asymptotes in the Ising regime. Although this effect occurs outside 
the finite-size regime, we may well hope that the so-called finite-range corrections applied in the previous section 
[Eq. (0)] can be used here as well. Indeed, these corrections are part of the universal scaling functions and although 
the amplitude bo = lim/j^oo limi^oo R 3 / 4 L 1 / 8 (\mL(K c )\) is a specific limiting value, the range-dependent correction 
factor does not depend on this limit. Especially the collapse obtained in Fig. g(b) makes it very tempting to apply a 
similar correction here. On the other hand, these corrections were calculated in the Ising regime, which we here are 
gradually leaving. In Fig. ^(b) we show the same data but now divided by the correction factor. Although a perfect 
collapse is not obtained, the asymptotes lie much closer together than without this correction. 

Also the critical amplitude of the Ising asymptote is known exactly. Indeed, by expanding Onsager's expression for 
the spontaneous magnetization |l7j , |l8|| , 
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sinh 4 (2 J/fc B T) 



1/8 

(9) 



G 



around the critical point J /k^T c — ^ arcsinh(l) we obtain for t < 

to = [4\/2arcsinh(l)(-<) + 0(t 2 )] 1/s m 1.22240995(-i) 1/8 . (10) 

For the nearest-neighbor Ising model R = R m = 1, so the fact that in Fig. ^ along the horizontal axis tR 2 is used 
instead of t and along the vertical axis (\m\)R instead of (|to|) does not affect the amplitude of the asymptote. 
However, the correction factor C[m] [denoting the factor between square brackets in Eq. @] must of course be taken 
into account. This correction factor describes the deviation of the critical amplitude do(R) from the leading scaling 
behavior in terms of a power series in R~ 2 (with coefficients that depend on In R) and it is not a priori clear whether 
C[m] converges for R = 1. It is certainly unlikely that a single term [the term proportional to 62 in (^]) vanishes] 
describes the deviation very well. No exact result for d (R = 1) = lirn^oo m^K^L 1 / 6 is known to us, but from a 
modest Monte Carlo simulation we found d a (R = 1) = 1.0092(4). On the other hand, from Eq. (|) with b Q = 1.466(2) 
and bi = —0.305(1) we find do(R = 1) = 1.018(4) which differs approximately two standard deviations from the 
numerical result. Recall that bo and b\ were obtained from a least-squares fit to the critical finite-size amplitudes for 
2 < R^ < 140. Nevertheless, the relative difference lies below the one-percent level, which cannot be distinguished 
in our graph. Therefore we have drawn the Ising asymptote with amplitude [4\/2 arcsinh(l)] 1 / 8 /(l — 61) in Fig. ^(b) 
and it indeed turns out to be a precise tangent to the crossover curve. 

As mentioned above, also the saturation effects can be described with mean-field theory. Namely, the magnetization 
follows from the well-known expression ||l9| , f2"of 

m = tanh(— ^-to] . (11) 



, T 

Rewriting this as to = (1 + t) arctanh(?7i) and solving for m one obtains below T c for small t 

to = V3M) 1 / 2 - \^{-tf' 2 - §- h ^{~tf/ 2 - y!V3M) 7 /2 + J|_ V3(-t)^ 2 + 0((-t)^ 2 ) . (12) 

The leading term shows the classical value (3 = 1/2 and the critical amplitude -\/3- To describe the saturation effects 
in Fig. the first three terms of this series suffice. Figure ^|(b) shows for the five largest ranges (i?^ = 140, 500, 
1000, 4000, 10000) the curves 



i\)R = Vsi-tR 2 ) 1 



/2 



1 - -^(-tR 2 ) - -±?—{-tR 2 ) 2 
5i? 2V ' 175i? 4V ' 



(13) 



For R m = 140 this expression does not precisely coincide with the numerical data, but for the remaining values the 
curves accurately describe the saturation effects. For these cases the interaction ranges are apparently large enough 
to suppress the critical fluctuations to a large extent. The lowest temperatures shown in the figure are T/T c = 0.52, 
0.60, 0.60, 0.68, and 0.50 for R 2 n = 140, 500, 1000, 4000, and 10000, respectively. Saturation effects become visible 
in Fig. ^| for t < —0.15, i.e. T/T c < 0.85. According to Eq. ( [l2|) the magnetization deviates here approximately five 
percent from the asymptote. Using Eq. (|l^) we can perform another operation on the numerical data. Namely, the 
influence of saturation effects in the mean-field model is described by the ratio between the full series expansion on 
the right-hand side of (U3) and its first term. As the mean-field expression constitutes an accurate description of the 
saturation effects for i?„ > 500, the factor between square brackets in Eq. (|l^) will give an accurate description of the 
relative saturation effects (i.e., the ratio of the saturated magnetization and the crossover curve) down to probably 
even lower interaction ranges. To illustrate this we have divided the data for R^ > 72 by the corresponding factor. 
The resulting graph [Fig. ^(c)] — in which also the data points in the finite-size regime have been omitted — shows that 
the data for these large ranges now nicely coincide on one curve, which is the actual crossover curve for the order 
parameter. 

The fact that for different interaction ranges the data (which overlap for considerable intervals of tR 2 ) coincide on 
one curve lends strong support to the hypothesis that the crossover curve is universal. As has already been noted in 
Ref. ]ll[, this contrasts with the conclusion drawn by Anisimov et al. || for the three-dimensional case. In Ref. Q 
it has been suggested that in the crossover region microscopic cutoff effects are not negligibly small compared to the 
finite correlation length £, which implies that the form of the crossover curve depends on the ratio between £ and the 
lattice spacing a. In our simulations we have not measured the correlation length directly, but we can still make a 
rough estimate from the data. Namely, at the locations marking the boundaries of the finite-size regime for different 
interaction ranges and system sizes in Fig. ^ the correlation length is approximately equal to the system size. From 
the magnetization densities for R^ > 72 we conclude that £ « 0.5/(— t), independent of the interaction range. The 
latter conclusion is in agreement with the above-mentioned renormalization prediction that the width of the finite-size 
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regime is to leading order independent of the interaction range. Thus, at a fixed value of the crossover variable tR 2 
the correlation lengths for different ranges have different values. However, the crossover curves coincide at fixed tR 2 
and hence are independent of the ratio £/a. 

Finally, we make some observations concerning the size of the crossover region. It is clear that it takes between 
two and three decades in the crossover variable to cross over from Ising-like to classical critical behavior. Thus, 
unless one studies systems with a rather large interaction range, one has to go to such a large temperature distance 
from T c to sufficiently decrease the correlation length compared to the interaction range that one has already left the 
critical region before observing classical critical behavior! The center of the crossover region lies in the neighborhood 
of |t-R 2 | = 1, consistent with a value for u of order unity. 



C. Magnetic susceptibility 

Unlike the order parameter, the magnetic susceptibility displays crossover upon approaching the critical point 
either from below or from above. We will discuss these two situations separately. In the ordered phase, T < T c , we 
approximate the magnetic susceptibility by the so-called connected susceptibility, 

& = xd K)-(IH) 2 . (14) 

In the two-dimensional Ising model with interaction range R this quantity will, close to the critical point, diverge as 
(— t)~ 7 / 4 i?~ 3 / 2 . Further below T c it will cross over to classical critical behavior, where x °c ( — I n a graph showing 
results for various ranges as a function of the crossover variable tR 2 a data collapse is obtained for x/ R 2 ■ However, 
just as for previous crossover curves, the data for small R will display an offset because of corrections to the leading 
dependence. To determine these deviations we first study the critical amplitude of the connected susceptibility, 
which was not considered in Ref. [^Cj| ; see Fig. [l^. The statistical uncertainty of this amplitude is notably larger than 
for (|m|) and (m 2 ) (cf., e.g., Fig. H), but one can still observe that the asymptotic regime is reached. In this figure we 
have also plotted the critical amplitude of the so-called scaled susceptibility k^Tx which was studied in, e.g., Ref. ||. 
Evidently, the latter amplitude shows a much stronger deviation from the leading range dependence, due to the fact 
that also T C (R) deviates from X" C MF (Fig. |l|). Thus, although both amplitudes have the same asymptotic behavior 
for large interaction ranges, it is much more difficult to extract this behavior from medium-range results for k^Tx 
than from the corresponding results for x- This may partially explain the difficulties experienced in Ref. ||. The 
deviations have been fitted to a correction factor of the form [1 + R~ 2 (v\ + V2 lni? 2 )], which we abbreviate as C\x\- 
Indeed, the fact that the finite-range corrections for x are so small allowed us to neglect them altogether in Ref. [nil , 
where only logarithmic scales have been employed. 

In Fig. pKa) we show the connected susceptibility, appropriately scaled with R and divided by the correction 
factor C[x]7as a function of the crossover variable. Just as for the magnetization density, deviations from the 
crossover curve are present even after the finite-range corrections have been applied. These effects are either caused 
by finite-size effects (close to T c ) or by systems that leave the critical region. In the latter case, saturation effects start 
to come into play. The finite-size effects are clearly recognizable in the rightmost part of the graph, where the curves 
start to follow horizontal lines. Once the temperature has been sufficiently decreased, the graphs start following an 
asymptote with slope —7/4, on which the data for various ranges quite accurately collapse. The amplitude of this 
asymptote is simply related to the exactly known amplitude A~ = 0.025537 . . . Q of the reduced susceptibility xo- 
This reduced susceptibility is defined as xo = ^b^x/^ 2 , where \x denotes the magnetic moment of a spin. This 
magnetic moment has also in our calculations implicitly been divided out. However, we should keep in mind that 
we have expressed all temperatures in terms of the mean- field critical temperature, i.e. T c = l/(zK c ), where z is the 
coordination number. For the nearest-neighbor model this yields an additional factor 4 and we thus expect a critical 
amplitude 2 arcsinh(l)j4~. In addition we have to take into account the finite-range correction factor which has been 
divided out. The question whether this factor is applicable for R = 1 has already been discussed in Sec. IV B| [below 



Eq. ([[(])]. Here, the difference between the deviation from the leading scaling behavior as predicted by C[x\ and the 
numerical result is approximately 3%, whereas the smallest differences that can be discerned on the logarithmic scale 
of Fig. |l^(a) are of the order of 5%. The asymptote with the above-mentioned amplitude divided by C[x] indeed 
lies tangential to the crossover curve, confirming our data. As the temperature is further decreased, the data for 
systems with small interaction ranges start to follow a line with a slope between that of the Ising and the mean-field 
asymptotes. This effect is caused by the fact that these systems have left the critical region. For sufficiently large 
interaction ranges, however, the curves coincide and have a slope that gradually decreases (in the absolute sense). 
Although the crossover curve at first varies more rapidly than for the magnetization density, it subsequently only 
slowly approaches the classical regime and the overall size of the crossover region is again between two and three 
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decades. Remarkably, the slope of the crossover curve passes even through the mean- field value —1 before settling at 
this value for sufficiently low temperatures. In other words, the derivative of the connected susceptibility appears to 
change nonmonotonically from its asymptotic Ising value —7/4 to its classical value —1. Several explanations may be 
considered for this behavior. Either it is an intrinsic effect of the crossover function or it might be explained from the 
fact that x is the difference between (m 2 ) and (|m|) 2 , which each separately are described by a monotonically varying 
curve. 

The saturation effects can — just as for the magnetization density — for large ranges be described with mean-field 
theory. In a mean-field model the magnetic susceptibility is given by 



1 



t + m 2 



(15) 



Using Eq. (|T|) we find for T < T c 



1 9 18 . . 18 , ,o 6714 , ., mf . ... 
*= ~t - 10 + 175 + 175 ( - £) + 67375 ^ + °^ > < (16) 

which exhibits the classical value for the susceptibility exponent, 7mf = 1 7 and the critical amplitude ^. Figure |TT|(a) 
shows the asymptote with this amplitude and one can observe that the crossover curve approaches this asymptote 
from below around tR 2 = -1. Also the mean-field curves @ are shown for R 2 m = 140, 500, 1000, 4000, and 10000 
and they accurately describe the numerical data. Thus, we have used the ratio between the series expansion ( |l6| ) 
and the asymptotic behavior 1 / ( — 2t) to remove the saturation effects in Fig. ^Jj(a). The resulting graph is shown in 
Fig. |llj(b), in which also the data points in the finite-size regime have been omitted in order to obtain a clear crossover 
curve. The nonmonotonical variation of the slope of this curve is clearly visible. 

In the disordered (symmetric) phase, we encounter a different situation. The susceptibility is now given by x' = 
L d (m 2 )/kBT. This is identical to the expression we ha ve use d for the finite-size crossover scaling, except that the 



temperature-dependent factor has been omitted in Sec. Ill C . Figure |l2j shows the critical finite-size amplitudes of 
both x' an d X — L d (m 2 ) as a function of the interaction range. We have fitted an expression of the form (Q) to the 
data for > 2. This expression describes the data well, except for the data point at R 2 n — 1, where the deviation 
is approximately 10%. Just as for the connected susceptibility, the finite-range corrections to the critical amplitude 
of x' are much smaller than for x- I n fact, they are so small that they can be completely omitted in the thermal 
crossover scaling, as illustrated in Fig. O. This graph shows x' ' /R 2 as a function of the crossover variable tR 2 for 
various interaction ranges and system sizes. Outside of the finite-size regime, the data follow the Ising asymptote 
with slope —7/4. The exactly known amplitude 2 arcsinh(l)j4 + , where A + = 0.96258 . . . ||lj, of this asymptote is 
accurately reproduced by the numerical data. For larger temperatures, the curves gradually approach an asymptote 
with the mean-field slope —1. However, some care has to be exercised when interpreting this behavior. Above T c , no 
saturation of the order parameter occurs and the system smoothly passes over to regular (noncritical) behavior. In 
this high-temperature region the susceptibility decreases proportional to 1/T. For small interaction ranges it is this 
behavior that one observes in the graph. Only for larger interaction ranges one actually observes classical critical 
behavior. The latter systems indeed reproduce the mean-field critical amplitude, which is equal to 1 [as follows from 



Eq. (15) with m = 0]. Note that, due to the absence of saturation effects, interaction ranges up to R^ = 1000 are 



amply sufficient to observe the full crossover region. 



V. EFFECTIVE EXPONENTS 

In several papers (see, e.g., Refs. |||§]) the slopes of the crossover functions are described by so-called effective 
exponents. These exponents can be defined as (3 c ff = dln(m)/dln \t\ = tdhi(m)/dt and 7 c ff = — c?lnx/c?ln |<| — 
—tdlnx/dt. In fact, this concept is familiar from the analysis of experimental data since a long time [ p2[ , but 
only a limited amount of theoretical work has addressed these issues. Of course, these exponents change from their 
Ising values to the classical values in the crossover region. However, the precise variation in the crossover region is 
unclarified and partially subject to debate. Although these exponents can be read off from the form of the crossover 
curves presented in the previous section, we consider it worthwhile to present separate graphs displaying /3 e ff, 7^, 
and 7^, where the superscripts denote the cases t < and t > 0, respectively. The additional advantage of these 
exponents is that they follow from data with the same range and hence are not affected by any range-dependent 
correction factors. 

As the graph in Fig. |2|(b) is particularly smooth, it is tempting to consider its derivative as well. As derived in 
Ref. jLO), {\m\)^/L oc (L/i? 2 ) 2 " 1-3 / 2 . This relation also holds in the mean-field regime, provided that one replaces the 







magnetic exponent by its starred counterpart yjj = 3d/4 = 3/2. The asterisk indicates that the exponent is modified 
due to the dangerous irrelevant variable mechanism, as explained in, e.g., Ref. p3| . Thus, while we can rewrite the 
above-mentioned relation in the Ising regime in terms of conventional critical exponents as (\m\)VL cx (L/ ' R 2 )-P/ u + 1 / 2 y 
this is not possible in the mean-field regime, since ^mf is not affected by the dangerous irrelevant variable mechanism. 
As an alternative we employ the specific heat exponent a; (|m|)\/i cx (L/i? 2 )~ 2 ^/( 2 ~ a ) +1 / 2 . The fact that the 
latter relation holds in the mean-field regime while the former does not is a direct manifestation of the violation of 
hyperscaling. Thus, we define [2/3/(2 — a)] e s = d\n({\m\)y/L)/d\n(L/R 2 ). This quantity is shown as a function 
of L/R 2 in Fig. [l4|. Although the error bars are considerable, the crossover from the Ising value 1/8 (for large values 
of L/R 2 ) to the classical value 1/2 (for small values of L/R 2 ) is clearly visible. 

Turning to thermal crossover, we display in Figs. [H| and ^ the exponents (3 e g , 7^, and 7^ as defined above. 
The effective magnetization exponent (3 c s increases monotonically from its Ising value 1/8 to the classical value 1/2. In 
particular do the data for different interaction ranges roughly fall onto the same curve, which supports the hypothesis 
that the crossover curve is universal. However, one observes that for systems with relatively small interaction ranges 
the effective exponent does not follow this curve. This effect, caused by saturation of the order parameter, can clearly 
lead to misleading results in experiments! In Fig. [l^ the nonmonotonical variation of 7^ between 7/4 and 1 is clearly 
visible. This may be considered as a manifestation of what Fisher Q calls an "underswing" . The occurrence of such 
a nonmonotonical crossover has been predicted by various renormalization calculations for the crossover from Ising 
to XY and Heisenberg critical behavior above T c , see, e.g., Refs. p4T-p6| and references therein. Furthermore, an 
exponent 7 c ff = 0.88(3) has been measured in the symmetric phase in micellar solutions |27|. Fisher || has suggested 



that an effective susceptibility exponent that takes a value 7 c ff < 1 in the crossover region might be a general feature 
of crossover from 3D Ising to classical critical behavior and noted that concrete calculations yielding such an effective 
exponent would be valuable. In Ref. [^J, a first-order e-expansion is quoted for the exponent crossover function, 

7eff = 1 + (7i - 7MF)£[ln(|t/G|)] , (17) 

where G is the crossover temperature or Ginzburg number and 

£(lny) = l/(l + ^/ 2 ) . (18) 

In our case, t/G is directly proportional to the crossover variable tR 2 . To describe the experimental results from 
Ref. [^7j, Fisher used the following extension of Eq. (|l8|), 

EQny) = (1 + P y e/2 )/[l + (p + l)y e/2 + qy £ ] . (19) 

Even though one may not expect such an expansion to converge for d — 2, we have drawn expression ([l7]) in Fig. |l6| , 
where we have taken the function E(hiy) from Eq. ([l9]), set e — 2 and adjusted p and q such that the curve 
constituted a reasonable description of the data. Clearly, no conclusions should be drawn from this curve, especially 
because Eq. ( p"9] ) has been proposed for the symmetric phase. In addition, for d — 2 the exponent e/2 is a very poor 
approximation for the exponent 9 = —y\/yt = 2, which is actually expected to appear in the function E{\ny). As 
follows from Fig. |l7], the behavior above T c is completely different. Here we have used expression (jl^) with Eq. ([l8]) to 
describe the data. Except for a shift along the horizontal axis (a proportionality constant in the Ginzburg number), 
no adjustable parameter is present and it is surprising how well the data agree with the theoretical prediction. 

Sometimes experiments have yielded effective exponents in disagreement with the known [ p8| universality classes, 
but still satisfying the scaling relations, such as 7 c ff + 2/3 c ff = 2 — a c s . Here a s denotes the effective exponent of the 
specific heat, which in our case is expected to be always (close to) zero, as both the classical and the 2D Ising value of 



a are equal to zero. This is also confirmed by the close resemblance between Figs. |14| a nd 15. Thus, it is interesting 
to note that this scaling relation is strongly violated in the present case: from Figs. |15| and |16| we can estimate that 
7 c ff + 2/3 e ff reaches a minimum of approximately 1.4 at tR 2 w — 1. 



VI. CONCLUSIONS 



In this paper we have presented numerical results for scaling functions describing the crossover from Ising-like to 
classical critical behavior in two-dimensional systems. While the general concepts describing this crossover have been 
developed many years ago, only a limited amount of progress has been made for a long time. In the present paper it 
is demonstrated, for the first time, that one can obtain accurate quantitative information on crossover scaling from 
computer simulations. The full crossover region was covered both for finite-size crossover and thermal crossover above 
and below T c . A data collapse has been obtained for all system sizes and interaction ranges, which supports the 
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hypothesis that these crossover functions are universal. Deviations from this curve are present but can be understood 
from finite-size and saturation effects. The results are in agreement with the previously derived renormalization 
scenario for these systems. 

Working in two dimensions offers the advantage that the exponents and the critical amplitudes are known exactly. 
More importantly, critical fluctuations are very large in two dimensions, which leads to critical behavior that strongly 
differs from classical behavior and hence to a clearly visible crossover between the two universality classes. We have 
shown that the magnetization density is described by a smooth crossover curve. The effective exponent, defined as the 
logarithmic derivative of this curve, increases monotonically from the Ising value to the classical value in two or three 
decades in the reduced temperature. On the other hand, the effective exponent for the susceptibility has a logarithmic 
derivative which varies monotonically above the Curie temperature and nonmonotonically below it. The occurrence 
of nonmonotonic behavior in the symmetric phase has been inferred from renormalization-group calculations in three 
dimensions and found long-standing interest. An extension of the present study to d — 3 is therefore highly desirable 
and has been planned for the near future. 
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APPENDIX: EXACT CALCULATION OF SOME FINITE-SIZE EFFECTS IN A MEAN-FIELD SYSTEM 

In Ref. Ejj ] the universal amplitude ratio Q has been calculated for a system in which all spins interact equally 
strongly, including the leading finite-size correction. It was shown that the relevant integrals can be expressed in 
terms of a quantity J/., which we here generalize to odd powers of the magnetization density, 



,1,,, ( »r^i-W = g) 1 V 



(Al) 



where N denotes the number of spins. In a similar fashion we can also calculate other finite-size effects to leading order 
in N. Expanding the terms in the partition function and replacing the sum over all possible states by an integral, we 
find the following expression for the absolute magnetization density 



(H) 



dm \di\ cxp j — — Nm? 



1 - — Nm 6 + -m 2 + (D{Nm 8 7 m 4 ) 
30 2 



(A2) 



in which Z denotes the partition function (except for a prefactor which has been divided out, cf. Eq. (31) in Ref. |L4j). 
Elementary algebra leads then to 



(\m\) - 12 1 / 4 ^^ 1 + O 



1 



N 3/4 



(A3) 



which for d = 2 yields expression (j^). Along the same lines one finds 



r(i) i 



12-^ 



m)AV2 



o 



N 



(A4) 



i.e. the susceptibility diverges as V N, and 



12 £(|)i 



o 



TV 3 / 2 



— + o . 

N V iV 3 / 2 



(A5) 
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FIG. 2. (a) Finite-size crossover curve for the absolute magnetization density multiplied by the square root of the system 
size, (b) The same graph but now the range-dependent corrections predicted by renormalization theory have been divided 
out. The correction factor abbreviated by C[m] stands for the factor between square brackets in Eq. (jfj]). A perfect collapse 
is obtained for all system sizes and interaction ranges. Both the exact mean-field limit (indicated by "MF") and the Ising 
asymptote with slope 3/8 are confirmed by the data. In this and all following figures the numbers in the key refer to values for 
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FIG. 3. Critical amplitude of (|m|) and the renormalization prediction fitted to it. This correction factor is used in Fig. ^(b). 
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FIG. 4. Finite-size crossover curve for the magnetic susceptibility divided by the system size. The range-dependent correction 
factor C[x] [the factor between square brackets in Eq. ([?])] has been divided out, as discussed in the text. Both the mean-field 
limit and the Ising asymptote (slope 3/4) are confirmed by the data. 
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FIG. 5. Finite-size crossover curve for the spin-spin correlation function multiplied by the system size. A range-dependent 
correction factor (abbreviated as C[g]) has been divided out, as discussed in the text. Both the mean-field limit and the Ising 
asymptote (slope 3/4) are confirmed by the data. 
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FIG. 6. Finite-size crossover curves for the amplitude ratio Q 



2 ) 2 /(m 4 ). Figure (a) shows the curve without any 



additional corrections, whereas in (b) a range-dependent correction factor C[Q] [see Eq. (Jm] has been divided out. For small 
values of the crossover variable L/R 2 the mean-field limit is reproduced and for large values of L/R 2 the Ising limit is approached. 
For a further discussion see the text. 
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FIG. 7. A detailed view of Fig. [](a) showing the deviations from the crossover curve for very small system sizes. 
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FIG. 8. The range-dependent correction factors C[m 2 ] — C[x], C[m 4 ], and C[Q] in (m 2 ), (m 4 ), and Q, respectively, as 
determined by least-squares fits to the critical amplitudes extracted from the Monte Carlo data. The line at height 1 is drawn 
for reference. One observes that C[Q] lies very close to, although not exactly at, unity. 
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FIG. 9. Thermal crossover for the absolute magnetization density for various ranges and system sizes. In figure (a) no 
additional correction terms have been used, whereas in (b) the factor C[m] has been divided out. In figure (c) the data for 
Rm > 72 have also been corrected for saturation effects and data points in the finite-size regime have been omitted. For an 
extensive discussion of the various features of these graphs the reader is referred to the text. 




FIG. 10. Critical amplitude for the connected susceptibility \ = L d ({m 2 ) — {\m\) 2 )/kBT as extracted from the thermodynamic 
limit of L~ 7 ^ 4 'xl(K c ). The dashed curve indicates the renormalization prediction fitted to the numerical data. Also the critical 
amplitude of the scaled susceptibility ksT c x is shown, which for small ranges deviates considerably stronger from the asymptotic 
behavior. 
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FIG. 12. Critical amplitude for the susceptibility \' ~ L d {m 2 )/kvT as extracted from the thermodynamic limit of 
Z/~ 7//4 Xl(^c)- The dashed curve indicates the renormalization prediction fitted to the numerical data. Also the critical ampli- 
tude of the scaled susceptibility UbTcX is shown, which for small ranges deviates considerably stronger from the asymptotic 
behavior. 




FIG. 13. Thermal crossover for the susceptibility x m tne symmetric phase, for various ranges and system sizes. No 
finite- range corrections have been applied. For a discussion see the text. 
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FIG. 14. The effective exponent [2/3/(2 — a)] e s as obtained from the finite-size crossover curve for (|m|)\/Z. 
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FIG. 15. The effective magnetization exponent /3 e ff describing the logarithmic derivative of the crossover function for the 
magnetization density. 



22 



1.75 



1.5 



1.25 



0.75 



2 
4 
10 
72 
140 
500 
1000 
4000 
10000 
8-expansion 



I ! ^ nL i - - ^ 



Ising 



-*§~*i 



MF 



" <i> 



-10- 



-10 7 



-10 



-10" 



-10" 



-10" 



tR z 



FIG. 16. The effective susceptibility exponent 7^ describing the logarithmic derivative of the crossover function for the 
connected susceptibility. The results on the left-hand side lie somewhat above the mean-field exponent due to saturation 
effects. 
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FIG. 17. The effective susceptibility exponent 7^ describing the logarithmic derivative of the crossover function for the 
susceptibility above T c . 
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p2 


2 


R 2 






500 
1000 
4000 
5000 
10000 


1580 
3148 
12580 
15704 
31416 


™*f » 251.770 
39 7 4 g 5 7 30 « 501.309 

1259568 w 2002.49 

98 392 7 6 59 « 2499 - 68 
"WTO 45 « 5000 - 34 


6.379(2) x 10" 4 
3.1904(6) x 10" 4 
7.9594(5) x 10" 5 
6.3746(3) x 10" 5 
3.18491(9) x 10~ 5 


6.3291139 x 10" 4 
3.1766201 x 10 -4 
7.9491256 x 10~ 5 
6.3678044 x 10~ 5 
3.1830914 x 10~ 5 



TABLE I. Some properties of the additional ranges used to span the full thermal crossover region. — 5000 has been 
included for completeness; it has only been used for the finite-size crossover scaling. The first three columns list the squared range 
of interaction R 2 ^, the corresponding number of neighbors z, and the squared effective range of interaction R 2 . Furthermore 
the critical coupling K c as calculated from Eq. (Q) and the mean- field approximation for the critical coupling K^ F = 1/z are 
shown. 
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